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Shape memory alloys inherit their macroscopic properties from their mesoscale microstructure originated 
from the martensitic phase transformation. In a cubic to orthorhombic transition, a single variant of marten- 
site can have a compatible (exact) interface with the austenite for some special lattice parameters in contrast to 
conventional austenite/twinned martensite interface with a transition layer. Experimentally, the phase compat- 
ibility results in a dramatic drop in thermal hysteresis and gives rise to very stable functional properties over 
cycling. Here, we investigate the microstructures observed in Ti 50 Ni5o-ccPd x alloys that undergo a cubic to 
orthorhombic martensitic transformation using a three dimensional phase field approach. We will show that the 
simulation results are in very good agreement with transmission electron microscopy observations. However, 
the understanding of the drop in thermal hysteresis requires the coupling of phase transformation with plastic 
activity. We will discuss this point within the framework of thermoelasticity, which is a generic feature of the 
martensitic transformation. 



I. INTRODUCTION 

Martensitic transformation in shape memory alloys (SMAs) is characterized by its macroscopic reversibility which gives rise 
to the shape memory effect (SME) or superelasticity d [2) . Macroscopic properties of SMAs strongly depend on their mi- 
crostructures originating from martensitic phase transformation (MT). MT is a first-order displacive transition consisting in a 
shear-dominated change of the underlying crystal lattice from a high symmetry phase, the austenite to a low symmetry phase, 
the martensite. The lattice mismatch between the different variants of martensite and the austenite results in long-range elastic 
interactions. Accommodation mechanisms such as twinning and tapering of martensite at the habit plane with austenite, mini- 
mize the strain energy. Despite these accommodation mechanisms, high levels of stresses are generated at the interface between 
the austenite and martensite. Several experimental studies have shown that dislocations were generated at the interface and that 
they build up during repetitive cycling through the phase transformation (3H3. Such dissipative, non-reversible mechanisms are 
responsible for the significant hysteresis observed in SMAs and the degradation of their functional properties over cycling which 
eventually leads to failure. The typical approach to tackle this problem is through thermomechanical treatments to strengthen 
the microstructure (strain hardening [8], nano- structuring OHO), precipitation hardening lfTTl[T2l ). 

Recently, however, a different approach has been investigated that simply consists in tuning the lattice parameters of some 
SMAs to obtain a geometrical fit between austenite and a single variant of martensite (phase compatibility). By doing so, a 
dramatic drop in hysteresis when phase compatibility is achieved has been observed H3][T4]|. Furthermore, in addition to a very 
small hysteresis, SMAs satisfying the compatibility condition have been shown to have very stable functional properties over 
cycling 1151 . 

So far, most of the SMAs studied that exhibit phase compatibility undergo a cubic to orthorhombic (c-o) martensitic transfor- 
mation giving rise to 6 variants of martensite described by their respective transformation stretch matrix . A matrix maps 
the austenite into one of the martensitic variants and its shape is dictated by the symmetry elements lost by the austenite. Ac- 
cording the geometrically non-linear theory of martensite (GNLTM) (0[T6]|) a compatible interface is possible when the middle 
eigenvalue A2 of the transformation matrix is equal to 1 . In most alloys this condition is not satisfied and therefore this is referred 
to as a highly non-generic condition in martensite. An example of a studied c-o system is TisoNiso-^Pd^ whose composition 
was systematically tuned to achieve geometric compatibility. In addition to a decrease in hysteresis, a dramatic change in the 
microstructure as one approaches A2 = 1 has been reported in reference [17]. The transmission electron microscopy (TEM) 
investigation of microstructures in these alloys shows the presence of large twinless martensite domains in contrast to internally 
twinned lamellar morphology commonly observed in martensites as illustrated in Figure [T] 

Different types of microstructures in alloys that undergo c-o phase transformation have successfully been studied in fT8l by 
means of the geometrically non-linear stress-free theories used for the microstructure optimization as proposed in references 
lfl~6l [T9l . It has been shown that one can obtain a very large variety of microstructure in c-o transitions. In a recent work, 
the formation of an exact austenite-martensite interface is observed by means of two- dimensional phase field method using 
geometrically linear elasticity (20). 

The goal of the present work is to carry out simulations of the phase transformation in SMAs using experimental data obtained 
for Ti5oNi5o- x ¥& x in order to reproduce experimentally observed microstructures. To that purpose, we use a three-dimensional 
phase field approach that incorporates the strain energy responsible for all specific features of the SMA such as transformation 
hysteresis and thermoelastic equilibrium. From there a discussion on hysteresis and thermoelasticity will follow. 



Figure 1 : TEM study displaying micro structure evolution with respect to Pd composition K2J\\22\l . (a) Internally twinned martensite plates are 
observed in Ti^Ni27Pd23 alloy with A 7^ 1, (b) in contrast, there are large twinless martensite domains in Ti^Ni^Pdn alloy with A 2 
very close to 1. 



II. PHASE FIELD METHOD 

The phase field method is a powerful tool for the investigation of microstructures at mesoscale. It has been successfully used 
in the study of a large variety of problems in solid-to-solid transformations l23ti25ll . The phase field approach consists in using 
some continuous fields of non-conserved and/or conserved parameters. These fields are functions of position and time and they 
are defined in the entire volume. They assume homogenous values within every phase domain but vary smoothly and contin- 
uously through narrow regions that correspond to phase boundaries. Martensitic phase transformations can be characterized 
completely by structural non-conserved long-range order parameters because these transformations are displacive solid state 
structural changes without atomic diffusion. The temporal evolution of the microstructure can be obtained by solving the usual 
dissipative time-dependant Ginzburg Landau (TDGL) equations: 

dt 6rj p (r,t)' 

where L is a relaxation coefficient, S denotes a functional derivative and r] p is the long-range order parameter associated to the 
martensitic variant p. 

In order to implement the evolution equations of the order parameters, one has to specify the free energy. The total free energy 
of a system in a volume V is written as follow 

hiatal = [ F GL dV + Elastic- (2) 
JV 

In the following, we describe each component of the free energy. 

The non-equilibrium chemical free energy of the stress-free crystal is given by a standard Ginzburg-Landau functional 

F GL {{ m }) = F L {{ m }) + § iv^l 2 , (3) 

i 

where Fl is an homogenous free energy and the second term on the rhs. is a heterogeneous (Ginzburg) term penalizing the 
spatial variations of the fields. The simplest fourth-order polynomial approximation of the Landau free energy density, invariant 
by permutation of martensitic variants, is given by 

m*}) = Aj2n> - Bj2v! + C(£ rtf. (4) 

iii 

This form may be used if the free energy is not invariant through the symmetry operations rji — » — rji. The coefficients A, B and 
C are positive temperature-dependent constants and can be fitted to reproduce the latent heat of transition, the equilibrium order 
parameters and the transition temperature. The gradient term prevents abrupt changes across the interfaces between different 
variants of martensitic phases and between the variants and the austenite. Thus, the order parameters change smoothly in a 
narrow region which means that the interfaces have a finite width in this formulation. Generally speaking the coefficient /? may 
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be chosen in such a way that the Ginzburg-Landau modeling reproduces some interfacial free energy density and/or interfacial 
widths. 

An arbitrary microstructure composed of different martensitic variants embedded in the austenite has a residual elastic energy 
(^elastic induced by lattice compatibility along the different interfaces. In the framework of linear elasticity, this compatibility- 
induced elastic energy is given by 

Elastic =\J C ijkl (ey (r) - e?.(r)) (e H (r) - e° kl (r))dV, (5) 

where Cijki are the elastic moduli BH . e^-(r) is the actual strain at a given coordinate r due to the deformation of the lattice 
and e?-(r) is the local eigenstrain. Thorough this paper, strain tensors and deformation gradients will be treated within the 
geometrically linear approximation. The stress-free tensor (p) of variant p is related to the corresponding transformation 
stretch matrix U (p) by 

4(p) = U(p)-I, (6) 

where I is the identity matrix. Within the linear geometry, compatibility between a martensitic variant and the austenite is 
possible when the middle eigenvalue of a stress-free strain tensor is 0. 

The actual strain (r) can be decomposed into homogeneous and inhomogeneous parts as 

€ij(r) = €ij + Ae^(r). (7) 

The homogenous strain e^ (r) describes the macroscopic shape of the system. Therefore, the inhomogenous strain Ae^ (r), 
which measures local deviation from the homogenous one, satisfies j v Ae^ (r)dV = 0, and may be related to an inhomogenous 
displacement field u(r) by 

_ l dudr) dujjr) 
A6 *'W- 2 1 g x . + dx . )■ 

Here, we consider the case where the crystal lattice remains coherent upon phase transformation i.e. the total displacement field 
is continuous. Furthermore, we assume that the elastic relaxation is much faster than the other relaxation mechanisms so that 
mechanical equilibrium is reached instantaneously. Therefore, we have: 

^elastic 



Sui(r) 



= 0. (8) 



These mechanical equilibrium equations can easily be solved in Fourier space for the displacement field. Inserting the solution 
back into Eq. [5] an explicit representation of the elastic energy follows: 

"V *V 1 /* * 6^q 

(^elastic = -^CijkiCijCki - VC ijk ieij (e kl ) + — C ijk i (e° kl e° kl ) - J n i cr 2 J (q)%/ c (n)cr^(q)n^ (9) 

where X(q) = y J v X(r)e - * q,r <ir is the Fourier transform of X, q is the wave vector, n = y^jj and (/(r)) = y J v /(r)dr. 

The tensor f^fc(n) is the Green tensor of the linear elasticity defined by Q~ k (n) = CijurijUi. The notations J* in Eq. prefers 
to the part q = that is excluded from the integration domain. The strain energy ^ depends on macroscopic deformation, 
which itself depends on the choice of boundary conditions. Thorough this paper, we adopt clamped boundary condition lij = 0, 
which is a convenient choice since in most of SMAs the macroscopic deformation is zero. 

In all generality, the local eigenstrain e^- (r) that enters into the elastic energy $ e iastic must be related to the long-range order 
parameters rj p (r) that enter into the Ginzburg-Landau functional . To the lowest order, and as here the total free energy is not 
invariant by the transformation r] p — )> we have: 



Using the above relation, the strain energy may be expressed explicitly in terms of order parameters and it is given by 

^ border 

^elastic 

z 

p,q=l 

The matrices B pq are defined as 



1 border r> i 

= 2 £ y y ^(q)^(q)^*(q)^3, dD 



B pq (q) = C ijk ie° ij (p)e° kl (q) - ^^(^^^(n)^^^)^, (12) 

where cr^-(p) is defined by a®- (p) = Cij k ie kl (p). The matrices B pq (q) characterize the long range elastic interactions between 
different phases. In equation flT] ), the notation X* is the complex conjugate of the quantity X. 
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III. NUMERICAL STUDY 



Our aim here is to perform dynamical simulations using the phase field model to understand the influence of the middle 
eigenvalue A2 on the microstructure. We perform simulations corresponding to identical physical systems in terms of thermody- 
namical potential and interfacial energy but that differ in terms of lattice parameters . We have chosen to use two TiNiPd alloys, 
namely TisoNiso-^Pd^ with x = 0.11 and x = 0.25, for which A2 = 1 and A2 = 1.018 respectively, to monitor the effect of 
the value of A2 on microstructure evolution (see reference l22l for lattice parameters). Additionally, we repeat the simulations 
using the lattice parameters of TiNiCu, satisfying the A2 = 1 condition but with an important volumetric strain (Tr(e^) w 0.02) 
associated with the martensitic transformation. The volumetric strains of TiNiPd alloys are much smaller i.e. Tr(e^ ) ~ 0.004. 



A. Simulation parameters 

The input parameters of the phase field method are the elastic constants, the transformation strains, the interfacial energy and 
free energy functional. To make sure that micro structural evolutions will depend only on the value of A2 for the three alloys in 
question, all simulation parameters are chosen to be same, except the transformation strains (e^.). 

The total energy density in a reduced form is written as 

F({V P }) = Fl + ^|V^| 2 + E elastlc , (13) 
where Fl = /? = -A12 , V = d Q V, E e i astic — ^ E ^ last ^ c with £ = C4 , 4£ ° . Here, fo is a typical energy density, eo is a typical 

JO JO^o ^44£q JO 

strain, do is a grid size. The coefficient £ is the ratio of typical strain energy to the driving force. 

In the simulations, we choose the following coefficients: A = 0.2, B = 3 and C = 2, fj = 0.5 and £ = 1, with these 
values, the order parameter rj p varies from zero in the austenite to 1 in the martensite. The typical strains are eo = 0.108 for 
Ti 50 Ni39Pdn, eo = 0.124 for Ti 50 Ni27Pd23 and eo = 0.0735 for TiNiCu. For each alloy, we use the same elastic moduli: 
Cn = 142GPa, C12 = 125GPa and C12 = 95GPa. The grid size of the simulation can be found using the relation between 
the experimental interfacial energy a and the dimensionless interfacial energy a 

a = af d . (14) 

One can introduce an interface between two variants and then allows the system to relax to its equilibrium without elasticity. 
The dimensionless interfacial energy a of such a configuration is given by 

<? = E{^- i? -/+Ef v i 7 ?pi 2 }' (is) 



where F re f is the free energy of a single variant reference state (without interfaces). Using Eq. 14 d is found to be 



do = ^ — 12- ( i6 ) 



Using the parametrization presented above, the grid size do is calculated to be 17nm, 21nm and 49nm for the alloys 
Ti 5 oNi3gPdii, Ti 50 Ni27Pd23 and TiNiCu, respectively. The required experimental interfacial energy for this calculation was 
taken from [ 26 ] due to lack of specific experimental data for the alloys investigated here. 

The kinetic equations Eq. [T]are integrated in parallel, using a semi-implicit Euler scheme in reciprocal space [ 60 1. Also, the 
choice of the numerical value of the stifness parameter $ ensures that the interfaces are large enough to avoid pinning on the 
numerical grid. 



Initial conditions 



The martensitic phase transformation has an athermal character. This means that thermal fluctuations do not play any major 
role. In general, the nucleation occurs around crystal defects [27 ] that decrease elastic energy barriers; whereas the forces due 
to thermal fluctuations are small compared to mechanical forces. Dislocations can be incorporated in the phase field model 
using the equivalence between the elastic fields created by a dislocation loop and those of a platelet inclusion [28-30]. A single 
dislocation was put in the middle of the system and all order parameters {7^} were set at zero at time t = 0. The dislocation loop 
lies at the yz plane and has a shape of a square with 3 grid points along every edge. Its normal vector is given by n = (1, 0, 0) 
and the Burgers vector is chosen to be b = (0, 0, 1). The simulation results are presented in the following. 
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Micro structural evolution in Ti^Ni^Pdn 



Figure [2] shows the microstructural evolution in T^oNisgPdn alloy. The figures correspond to the middle plane perpendicular 
to the z direction of a 256 x 256 x 256 simulation box. The temperature was held constant during the entire evolution. The 
austenite (black) transforms into variants of martensite around the dislocation loop. All six variants appear simultaneously 
(variants 1, 2, 3, 4, 5 and 6 are shown by red, green, blue, white, yellow, light blue, respectively). It is important to point out 
that we do not observe the growth of any typical poly twinned plates formed by pair of variants with specific twin ratios. Rather, 
we observe a chaotic microstructure formed by all twin-related variants with inequivalent volume fractions. For example, it can 
be seen from figs. ^ fo)||(c) that there is a rather large region formed by variant 4 (white), while other variants occupy smaller 
volume fractions. The simulated microstructure is qualitatively similar to the ones observed experimentally, as seen in Fig. [3] 
where we present a comparison between simulation and observation. It is clear that variants with large or small volume fractions 
are present in the crystal as well as in the simulated microstructure. 

The origin of this chaotic microstructure is essentially due to the fact that the A2 = 1 condition allows for the growth of 
monovariant domain within the austenite without the need of a twinning mechanism to relax stress accumulation along the 
martensite/austenite interface. The transformation process may be qualitatively described as follows. 

As mentioned previously, the nucleation occurs on single defects, here a simple dislocation loop. The first embryo is inho- 
mogenous, but does not display any regular stacking of twins. Its morphology and initial growth are not controlled by compat- 
ibility relation between martensite and austenite, nor by the twin relations between martensitic variants, but by the stress-field 



of the pre-existing defects (see Fig. |2(a)| ). At this stage, the interfaces between the embryo and the austenite do not necessarily 
fulfill the compatibility relation. Any component of the embryo interface which is non-compatible will favor, by a local twinning 
process, the subsequent growth of other variants. If, however, it happens that at some point of the formation process, a variant 



displays a nearly compatible interface with the austenite, as variant 4 (white) in Fig. |2(b)| it may adopt a "pre-needle" shape 
composed of two nearly parallel and stress-free interfaces that terminate on a rounded tip. In front of the tip, the stress in the 
austenite will be high and, therefore, the transformation driving force much higher than along the two nearly compatible sides of 
the needle. This will favor the propagation of the tip and therefore, the rapid growth of the needle-shaped martensite, as seen in 
Figs. 2(a)|2(d) where we observe that variant 4 rapidly reaches the system boundaries, and therefore, occupies a larger region 
than the other variants in the final microstructure. 

In order to understand the absence of laminates in the microstructure, we checked that, when the A2 = 1 condition is fulfilled, 
the residual elastic energy associated to a poly twinned domain formed by two different variants is always non-zero if the volume 
fraction of the variants are finite (see Appendix). This analysis shows that a twinned-martensite/austenite interface will be 
energetically non favorable when \ 2 = 1. Hence, a polytwinned domain cannot grow within the austenite. 

The three dimensional final state of this simulation is shown in Fig. ^^a) Another microstructure of interest identified in this 
figure is the triangular structure (encircled region in Fig. |^a)| ) formed by three variants. The magnification of this structure is 



shown in Fig. f a) Similarly, one can see another self-accommodation unit in Fig. ^ a) formed by variants 4-1-3 (white-red 



blue). This kind of microstructure was also observed in the experiments when compatibility condition is fulfilled, as seen in Fig. 

HI 

In brief, we understand that the martensitic transformation may proceed by using simultaneously (along different portions of 
the austenite-martensite interface) or alternatively two mechanisms: a "local twinning" mechanism as well as the "compatible" 
one. The additional freedom afforded by the later is at the origin of the chaotic morphology developed by the martensite, at the 
expense of the more regular polytwinned domains that would be observed if the twinning mechanism was operating alone. 

















(a)t = 23 x 10 3 



(b)t = 29 x 10 3 



(c)t = 30 x 10 3 




40 x 10 3 



Figure 2: Time evolution of the microstructure for Ti^Ni^Pdw alloy, for which A2 — 1. The simulation box size is 256 x 256 x 256. 
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Figure 3: Comparison of the simulated micro structure with experimental observations 07V . In both cases, A2 = 1 condition is satisfied. 




(a)t = 55 x 10 3 



(b)t = 168 x 10 3 



Figure 4: Comparison of the final states of the alloys with different A2. (a) Final state of Ti^Ni^Pdw (X2 = I), and (b) 7150^27^23 




Figure 5: (a) Magnification of the encircled region in Fig. \4j\a)\ It shows a typical self-accommodation unit formed by variants 6-3-2. (b) TEM 
image displaying the micro structure ofTi^Ni^Pdn alloy. It exhibits a triangular self-accommodation unit similar to the simulation results. 



Micro structural evolution in 7/50^/27^23 



Figure [6] shows the microstructural evolution in Ti 50 Ni27Pd23 alloy. We use the same initial conditions as in previous simu- 
lation. As above, the first embryo is inhomogenous and its initial growth is controlled by the stress-field of the dislocation (see 
Fig 6(a)| ). After this initial growth, poly twinned domains compatible with the austenite appear. For instance, variant 1 (red)- 




Figure 6: Time evolution of the micro structure for Ti5oNi27Pd23 alloy, for which A2 7^ 1. The simulation box size is 256 x 256 x 256. 




Figure 7: Time evolution of the micro structure for TiNiCu alloy, for which A 2 = 1 but the volume change associated with the transformation 
is non-null. The simulation box size is 256 x 256 x 256. 



variant 5 (light blue), variant 6 (yellow)-variant 2 (green) and variant 6 (yellow)-variant 3 (blue) form laminates that grow into 
the austenite as shown in Figs. ^ fo)||(d) The comparison of Fig. 2(d) and Fig. 4(a) that correspond to a situation with A2 = 1, 
with Fig. 6(d) and Fig. 4(b) , for which A2 = 1.018, shows that the micro structure evolution is very sensitive to the value of 



the middle eigenvalue A2 : a slight deviation from compatibility generates a micro structure which is significantly more regular 
than the one observed if A2 = 1. Finally, in Fig. ^^b) we show the three dimensional final state of the microstructure, where we 
observe several twinned laminates in contrast to the final state of Ti 50 Ni39Pdn alloy. 



Micro structural evolution in TiNiCu 



The microstructural evolution in TiNiCu alloy satisfying the condition A2 = 1 but with an important volumetric strain in the 
martensitic phase subjected to same initial conditions as in previous examples is shown in Fig. [7] However, the temperature was 
decreased continuously. The reason for this is that the presence of the large volume effect causes the transformation to reach 
its thermoelastic-equilibrium. As previously, in the early stage, the complex multi-variant structure grows into the austenite 
around the defect. But even at late times, there is retained austenite in the system. We observe that variants form two nearly 
parallel compatible interface with the austenite that terminate on a tip. As explained previously, in front of the tip, the stress in the 
austenite is high, and thus, it results in higher transformation driving force. Therefore, mono-variant domains grow perpendicular 
to the habit plane normal direction. When a variant tip reaches the system boundaries and/or meets another martensitic domain, 
the further growth in front of the tip is not possible. At this stage, we observe that variants grow in the direction of the habit 
plane normal as shown in Figs. H[c)|(d)| in the encircled region. There is no self- accommodation by twinning as observed 
experimentally (32). For this alloy, microstructural evolution is affected by both the large volume effect and the value of A2. 



IV. THERMOELASTICITY 



The phase field simulations presented above reproduce microstructures in good agreement with experimentally observed 
microstructures. However, when we calculated the volume fraction of martensite with respect to the temperature, we did not 
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find a significant influence of the A2 value on the transformation hysteresis. In this section, we discuss this contradiction by 
considering the thermodynamics of thermoelastic martensitic transformations. 

Thermoelasticity is a generic feature of the martensitic transition in SMAs. Basically, it manifests itself by the fact that the 
martensitic transition spreads over a finite temperature range, generally with a small thermal hysteresis. The forward transition 
(from austenite to martensite) starts at a temperature M s and generally finishes at a lower temperature Mf . The reverse trans- 
formation on heating starts at a temperature A s and in general comes to completion at a higher temperature A f with Mf < A s 
and M s < Af. The width of the hys teresi s may be quantitatively characterized by | M s — Af \ and | Mf — A s | and its spread by 
\M f - M s \ and \A S - A f \ (see Fig.[9(b)|). 

In principle, in the context of equilibrium thermodynamics, the martensitic transition should take place at the temperature 
To for which the stress-free (or chemical) Gibbs energies of the austenite and martensite, given by Gm = Hm — TSm and 
Ga = Ha — TSa (i.e., the free enthalpies), are equal. 

If the transition did take place precisely at the equilibrium temperature To, it would be completly reversible, with no hysteresis 
and accompanied by a latent heat AQ (measured by calorimetry) equal to the difference between the austenite and the martensite 
enthalpies: AQ = H A (T ) - H M (T ). 

In reality, the situation is more complex and, as mentioned above, the forward and reverse transformations spread over a finite 
temperature range and display a (small) hysteresis. A large quantity of experimental and theoretical studies have been devoted 
to understand the origin of those two phenomena (see (33), |[34l and ifSTll for reviews). 

However, there is still often some confusion concerning the correct comprehension of these features, even though it is clearly 
understood that the origin of thermoelasticity and of the finite width of the hysteresis are due to stored elastic energy and 
frictional work. 

The underlying physical mechanisms responsible for the spreading of the forward and reverse transitions and of the hysteresis 
are qualitatively described in the following. We in particular clearly explain that thermoelasticity by itself is a thermodynamically 
reversible process (in other words, the process of storing elastic energy is not irreversible), whereas the presence of local free 
energy barriers (such as those due to friction on local external defects or to the competitive growth of different martensitic 
variants) leads to irreversibility and intermittent dynamics. This analysis is illustrated through the discussion of two different 
simple situations that involve the interplay between stored elastic energy and a local constraint (the model of Olson and Cohen) 
or fixed boundary conditions (a clamped system). 



The martensitic transition in a perfect single crystal 



Due to the large strain associated to the martensitic transition, nucleation of a martensitic domain requires accommodation by 
twinning. Therefore, the first embryos are essentially modulated or simply sheared as onset for the poly twinned character of the 
final martensite domain, and cannot nucleate homogeneously through thermal fluctuations [ 35 , 36 ]. Nucleation is heterogeneous 
on pre-existing defects that generate local favorable strains. Note that, to a certain extent, these embryos exist above the equilib- 
rium temperature T and can be responsible for the softness of the effective shear elastic constant often observed experimentally 
well above To. From this point of view, the "soft mode" phenomenon often referred to as pre-martensitic and cited to be a 
characteristic of the martensitic transition, is not an intrinsic softening of the austenite phase, but rather a collective and effective 
response of pre-existing martensitic-like embryos stabilized by static defects, such as dislocations and/or grain boundaries or 
even pinned defects (35j [36j [42) 

Inside the martensitic embryos accommodation cannot be perfect, because the boundary of any finite domain cannot consist 
exclusively of habit planes. Therefore, elastic screening by twinning is only partial and some unrelaxed elastic energy is stored 
in the resulting martensitic domain and in its surrounding. Below To, when Gm(T), the Gibbs free energy of martensite is lower 
than Ga(T), the Gibbs free energy of the austenite, part of the stress-free driving force AG(T) = Gm(T)—Ga(T) is consumed 
by the stored elastic energy. A simple model (see below) shows that, for sufficiently small undercooling, a modulated embryo 
that nucleates on a defect reaches an equilibrium state with a finite volume. When temperature is sufficiently lowered, this 
embryo becomes unstable and undergoes a continuous growth at a fixed temperature. In brief, the martensitic start temperature 
M s is lower than T , and corresponds to the first embryo that becomes unstable. M s is not intrinsic and it depends on the 
nature of the pre-existing defects, and may change during cycling (because the martensitic transition itself creates defects). The 
growth of this isolated embryo is irreversible, because it occurs at a temperature when the martensite and the austenite are not 
in thermodynamic equilibrium. Therefore, it is accompanied by an increase of the total entropy (entropy of specimen plus the 
entropy of the surroundings). If the embryo does not meet any obstacle, the growth is continuous and, in an ideally perfect 
single crystal, the martensite growth will reach completion at Mf = M s : no thermoelasticity would be observed, even though 
the growth is continuously accompanied by a stored elastic energy that consumes part of the available stress -free driving force. 
Since the martensitic phase is more stable than the austenite for temperatures smaller than To, the reverse transformation on 
heating will start at a temperature A s strictly above M s . 

However, A s may be below or above T depending on the amount of unrelaxed elastic energy stored in the microstructure 
when completion is reached and the temperature reversed. More precisely, A s is the smallest temperature at which, on heating, 
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some martensitic domain becomes unstable with respect to the austenite. If the specimen under consideration is a prefect 
single crystal and if the forward transition involved only one martensitic plate, it is likely that, at completion, the stored elastic 
energy was essentially reduced to zero because an almost perfect self-accommodation may be reached when there is no residual 
austenite and if the external boundaries are free. As a result, the reverse transformation cannot take place before the stress-free 
bulk austenite become more stable than the stress-free bulk martensite. This implies A s > T . Also, as soon as the reverse 
transition starts, it will go to completion at fixed temperature, because all the components of the total Gibbs free energy operate 
in the same direction (the free energy difference, the stored elastic energy and the interfacial energy are all positive). Therefore, 
we have A s = Af. 

In summary, in an ideally perfect single crystal and if only one martensitic plate is involved in the transformations, we have 
M s = Mf < To < A s = Af. The forward and reverse transitions involve stored elastic energy but there is no thermoelasticity. 
The transitions are irreversible. 



Thermoelasticity with internal constraints 

The situation is very different if the transitions involve some constraints that may impede the growth or the shrinkage of the 
martensitic domains, such as grain boundaries or pre-existing precipitates of a secondary phase. These constraints may be due 
also to the martensitic transition itself if many differently oriented martensitic plates are involved, or if some specific boundary 
conditions are imposed, such as a fixed volume in the presence of a dilatational component in the eigenstrain of the martensitic 
phase. 

As we explain below using two very simple generic models, these morphological constraints may lead to a situation where 
a martensitic plate (if the constraints are local, such as due to grain boundaries...) or the martensitic volume fraction (if the 
constraints are global, such as fixed boundary conditions...) reach a stable equilibrium state where the driving force due to the 
stress-free Gibbs free energy difference between the martensite and the austenite is exactly equilibrated by the driving force due 
to the stored elastic energy. In such state, the volume of the martensitic plate (or simply the martensitic volume fraction in the 
second situation) is proportional to \Ag(T)\, where Ag(T) = g m (T) — g a (T) is the difference between the stress-free Gibbs 
free energy densities of the martensite and the austenite. In other words, in all cases, the volume fraction of martensite varies 
continuously with the quantity of undercooling: this is precisely the definition of thermoelasticity. An important feature of this 
state is that the stored elastic energy consumes only a part of the available stress-free Gibbs free energy (exactly half in the two 
generic models presented below). If no other phenomena come into play (such as local free energy barriers due to the friction 
of the moving interfaces or due to the competition between differently oriented martensitic plates...), the martensitic volume 
fraction varies continuously with temperature: the transition is reversible, even though stored elastic energy is involved. In other 
words, in its simplest expression, thermoelasticity by itself is a reversible phenomenon. 

Irreversibility and intermittent dynamics 

However, the situation described above is oversimplified. In real systems, the moving interfaces encounter local free energy 
barriers of different nature. First, there is the presence of defects such as dislocations already present before the transformation or 
created during the process itself, either by partial plastic accommodation of the unrelaxed elastic energy along the habit planes or 
by the twinning mechanism itself. Second, there is the competitive growth between different martensitic plates. In all generality, 
the moving interfaces must overcome these defects: this is the origin of the pinning phenomena and the associated frictional 
work . In term of the free energy landscape, these defects are associated to local energy barriers that cannot be overcome by 
thermal fluctuations. Therefore, for each event a finite amount of undercooling must be devoted to the decrease of the associated 
energy barrier until it disappears. As soon as the barrier vanishes, the system becomes unstable and moves, at fixed temperature, 
until it meets the next barrier and becomes again pinned. From the thermodynamics point of view, the evolution at decreasing 
temperature between two consecutive jumps, that takes place at Ti and T 2 < Ti, is reversible l43ll (the total entropy, i.e., the 
entropy of "universe", stays constant) and the release of latent heat AQ during this process (i.e., measured by a calorimeter 
during a calorimetry experiment) is equal to the negative of the variation of enthalpy: AQ(T\ — » T 2 ) = H(T\) — H(To) {44\ 
(here and below, it is meant that the enthalpy H, and therefore the free enthalpy G = H — TS, includes the elastic energy). 
When temperature reaches T 2 , the barrier collapses and the system jumps at fixed temperature to the next metastable state. This 
corresponding evolution is irreversible, because the ending and initial states do not have the same Gibbs free energy: the entropy 
of universe increases by an amount AStot given by AS to t(T 2 ) = - AG ^ 2 \ where AG(T 2 ) is the variation of the total Gibbs free 

energy. The jump at T 2 is of course also accompanied by a release of latent heat AQ(T 2 ) given by AQ(T 2 ) = #(T 2 ) — #(T 2 ) 
(here and below, it is meant that the enthalpy H and therefore the free energy enthalpy G = H — TS, includes the elastic energy). 
Part of this latent heat is released irreversibly. This irreversible latent heat is precisely equal to the negative of the variation of the 
total Gibbs free energy during the jump: AQ irr = — AG(T 2 ). If the system was coupled to a mechanical device, the irreversible 
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component of the released latent heat would correspond to "lost of useful work" due to irreversibility, i.e. to the part of the Gibbs 
free energy variation that is not converted into "useful work" (work other than PV). 

In summary, thermoelasticity during the forward and reverse transformations can be observed if and only if these processes 
are accompanied by a finite unrelaxed stored elastic energy and if some internal (local) or morphological (global) constraints 
impede the growth or the shrinkage of the martensitic domains. Without local free energy barriers, such as those due to defects 
(usual "frictional work") or to the competing growth of different martensitic plates, the thermoelastic evolution would be thermo- 
dynamically reversible. In real systems, these barriers exist and the evolution consists of a series of small reversible quasi-static 
paths between consecutive irreversible jumps at fixed temperatures: thermoelasticity is then globally irreversible. 

Such intermittent behavior has been observed in many slowly driven physical systems operating in nonequilibrium steady 
regimes associated, for instance, with magnetism, superconductivity, porous flow, friction, plasticity, fracture, earthquakes (50l 
[52l[53j[55l. These systems evolve toward a critical state represented by the ensemble of metastable states through which the 
system passes via an avalanche-like dynamics whose amplitudes and durations follow power laws (see Ref. l54l for a review). 
Similarly, acoustic emission experiments as well as calorimetry measurements in temperature-driven athermal martensitic phase 
transformations have shown that the martensitic transition proceeds through avalanches (see for example (48] [57l |62) , and also 
Ref. [ 63 1 for an interesting discussion on the influence of the driving rate with respect to the role of an external disorder). 



Model of Olson and Cohen 

The origin of thermoelastic behavior in martensites was first examined by Olson and Cohen l37l using a very simple ther- 
modynamical model. In this model, the free energy change (AG) of a growing ellipsoidal poly-twinned martensitic particle of 
radius r and thickness c was expressed by 

AG = ^-7rr 2 cAg ch + 27rrc 2 fiel + 2nr 2 a, (17) 

where Ag c h(T) = g m (T) — g a (T) is the difference between the Gibbs free energy densities of stress-free bulk martensite and 
austenite, /i is a shear modulus, eo is a typical stress-free strain and a is the interfacial energy density. The first term at the r.h.s. 
of equation ( [T7] ) corresponds to the stress-free driving force. It is supposed here E51 that the martensitic particle is polytwinned, 
in order to accommodate its bulk elastic energy. However, as already mentioned previously, this accommodation mechanism 
cannot be perfect. First, the boundary at the periphery of the ellipsoidal particle does not correspond to a habit plane. Therefore, 
this ribbon-like boundary generates a strain field which, if c « r, is equivalent to the strain field generated by a dislocation. The 
corresponding unrelaxed elastic energy l46l is given by the second term of equation §FI\ . Second, the width of the twins along 
the habit planes is not infinitesimal and, consequently, there is a stray stress generated by the alternating domains. This stress is 
localized near the habit planes and thus the associated unrelaxed elastic energy is proportional to their areas. This contribution 
is included into the interfacial energy density a, which in turn governs the last term of equation ( [T7] ). 

The contour plot of the free energy density is shown in Fig. [8] It is seen that the stable and unstable regions are separated by 

the dashed line passing through the saddle point with temperature dependent coordinates r s = A ^ e ^ 2 and c s = — Ag 2 ^ T ^ • At 
sufficiently small undercooling, the polytwinned domain is metastable, with finite radius and thickness. However, for a critical 
undercooling, the unstable region will move to the left due to temperature dependence. Therefore, the precipitate coordinates 
r ini and c ini will be in the region (right side of the dashed curve shown in Fig. [8]) where the precipitate is unstable: the growth 
of a particle cannot be stopped in the absence of any morphological constraints such as grain boundaries, defects, fixed external 
boundaries or similarly other transformed domains of different orientations. When radius and thickness of martensitic particle 
will be in the unstable region after the required driving force is supplied, the martensite volume fraction will reach 100% at 
the transition start temperature without increasing the driving force. On the other hand, if the radial growth is stopped by any 
of these morphological constraints at a given radius r = r max , the transformational force in the thickening direction leads to 
an energy balance condition between the stored strain energy and chemical driving force given by: This lead to an equilibrium 
thickness given by: 



- &9ch(T ) n 

3/i6Q 



^eq(T^ — 2 ^max- (1$) 



When this state is reached, the stored elastic energy AG e i = 2i\r rnax (? eq ii^ ) is exactly opposite to half of the available stress-free 
Gibbs energy AG ch = ^r^^c^Ag^: 

AG el (T) = -^^. (19) 



As seen in equation (18), the thickness of the martensitic particle beyond the critical undercooling varies continuously with 
temperature. This is the definition of the thermoelasticity. This simple model shows that the required condition for thermoelastic 
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Figure 8: Contour plot of the free energy density with respect to radius r and thickness c of the martensitic particle. The saddle point 



coordinates are given by: r sa ddie/ro = § and c sa ddie / 'c = f , where r = 4A 9 ^ ( ^ )2 and 



co 



3a 



2A 9ch (T)- 





Figure 9: Hysteresis of the growth of a single martensitic particle (a) in the absence of any kind of friction, (b) with friction, r is the volume 
fraction ofmartensite and T is the temperature. 



equilibrium is partial self- accommodation that leads to the stored elastic energy and internal constraints. The corresponding 
hysteresis of the growth of a single martensitic particle is shown in Fig. [9] The transition starts at a martensite start temperature 
M s below the transition temperature To defined by Ag c h = 0. This undercooling is required in order to overcome the energetic 
barrier due the interfacial energy and to the stored elastic energy, and in order to reach the saddle point. Afterwards, the 
martensitic particle grows along the segment (a-b) until its radial growth is stopped by a morphological constraint. At point b 
the particle reaches its thermoelastic equilibrium and along the segment (b-c) its growth is dictated by equation ( 18 ), i.e., the 
transformation is frozen if the temperature is not decreased. On the other hand, when the system is heated, the martensitic 
particle retracts along the segment (c-d). It is easy to argue that, on heating, point d is beyond point b, i.e., Af > M s , and that 
Af is lower then the equilibrium temperature To . Indeed, on heating, the martensitic plate whose radius r max is imposed by 
some internal constraint, becomes unstable when the radius of the saddle point, which is temperature dependent, becomes equal 
to r max . Obviously, the relation r ini < r max < oo implies M s < Af < T . Also, as the particle state at point d corresponds to 
an unstable state, the shrinkage at point d cannot be stopped. This scenario leads to the behavior displayed in Fig. |^a)| 

The hysteresis described within this model (see Fig. |^[a)| ) does not take into account any frictional stress for interfacial 
motion and it is reversible for temperatures below the critical undercooling, i.e., along segment (c-d) in Fig. ^^a) However, the 
frictional forces are always present in a system (see the discussion above) and create local and small energetic barriers that must 
be overcome in order for the interfaces to move. This mechanism introduces metastable states in the system and therefore, the 
hysteresis shown in Fig. |£|[a)| becomes as shown in Fig. |^b)| which is constituted by small bursts between these (meta)stable 



states, each of them at the thermoelastic equilibrium. 



Thermo elasticity for a non-volume preserving transformation and fixed boundary conditions 



In this subsection, we will show using a simple qualitative model that a non-volume preserving transformation together with 
fixed boundary conditions (zero macroscopic deformation or hard device boundary conditions) will lead to a thermoelastic 
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equilibrium in the absence any other morphological constraints. 

To start with, we consider a system formed by the austenite and polytwinned martensitic domains with a non-zero dilatational 
part. We assume that the strain energy due to the shear components of the strain tensor is completely accommodated. But, 
of course, the dilatational part cannot be accommodated. The stress-free strains of martensite and the austenite are eo and 
0, respectively. However, due to fixed boundary conditions, the actual strains e\ and 62 experienced by the martensite and the 
austenite will differ from the stress-free strains. We look here for a very qualitative picture. Therefore, as far as the elastic energy 
is concerned, we suppose that the strains e\ and 62 are homogenous. The total elastic energy AE strain is then approximately 
given by: 



KV =V {ei ~ eo) + y' 2 ' 



(20) 



where K is a bulk modulus, V is the total volume of the system, and V\ and V2 the volume of the martensite and the austenite, 
respectively. For a macroscopic deformation imposed to zero, one can write 



eiVi +e 2 V 2 = 0. 



(21) 



The above equation can also be rewritten in terms of volume fraction of each phase as 

em + e 2 r 2 = 0; 



(22) 



where t\ and r 2 are the volume fractions of the martensite and the austenite, respectively. The strain energy in terms of the 
volume fractions becomes 



AE S 



KV 



- e ) H e l5 

T 2 



(23) 



where we used e2 = — The minimization of the above strain energy with respect to ei assuming that the strain energy 



AE, 




Figure 10: Total free energy as a function of martensitic volume fraction for a non- volume preserving transformation under clamped boundary 
conditions (zero macroscopic deformation). 



quasi-statically reaches its mechanical equilibrium expressed by dAE *train _ q ? i eac [ s to 



dei 

ei = r 2 e . (24) 



By substituting the above equation into equation (23 ), one finds 



^^strain 2 2 /oc\ 

kv =Tl€ °- (25) 

The total free energy change AE to t including the stress-free Gibbs energy is then given by 

AE tot = (jle\K - n \Ag ch \)V. (26) 
This free energy change (see FigfTOb has a minimum for 



eq _ \^9ch 

fl 2Ke? 



(27) 



o 
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The above equation expresses the fact that for a given undercooling |A^ C ^|, the martensite volume fraction will reach the value 
of T^ q . Afterwards, the volume fraction of martensite can only further increase (continuously) by supplying extra driving force. 
This is the definition of the thermoelastic equilibrium. The value of the total free energy in the thermoelastic equilibrium is 
found to be 



AEZ = "^if v > (28) 



and again, as in the previous situation, the stored elastic energy consumes exactly half of the available chemical free energy: 

Strain = ~\t? Ag ch V = -\&E% (29) 

where A£^ is the chemical free energy change associated to the formation of the martensitic volume fraction 



For any further quasistatic temperature decrease 5T, Eq. 29 indicates that the increase in stored elastic energy SAE^ q 



strain 



always consumes exactly half the variation SAE^ of the chemical free energy made available by the increase of the martensitic 
volume fraction. In other words, only part of the available chemical free energy change is stored in the form of unrelaxed strain 
energy. However, eventhough the "internal work" stored in the transforming specimen is less than the available chemical free 
energy change, the thermoelastic transformation is thermodynamically reversible, simply because it consists in a succession of 
quasistatic equilibrium states. Because of the variational character of the total free energy at thermodynamic equilibrium, these 
states, correctly defined by the successive minimizations with respect to the strain and to the martensitic volume fraction in 
the present simplified model (or through the energy balance expressed in Eq. ?? in the previous example), correspond to true 
equilibrium states, and not to states where the equilibrium would only be established along the interfaces between the martensitic 
domains and the austenite, as suggested in (33). Also, further growth of the martensite is continuous and do not necessitate any 
finite increase of undercooling. 

In conclusion, this simple model shows that in the absence of any other morphological constraints, the volume change asso- 
ciated with the martensitic transformation together with the macroscopic deformation imposed to zero will cause the system to 
reach a reversible thermoelastic equilibrium. 



V. CONCLUSION 



In this study, microstructural evolutions of several Ti-Ni-Pd and Ti-Ni-Cu ternary martensitic alloys have been investigated. 
For specific concentrations, these alloys satisfy a very non-generic condition such that the middle eigenvalue A2 of their transfor- 
mation matrices is equal to 1. This condition allows the formation of compatible austenite/single variant of martensite interfaces. 
First, we realized a phase-field study on this subject and obtained different microstructural evolutions that depend on the value 
of A2. The microstructures obtained were in good agreement with the experiments, i.e., twinless domains are observed when 
A2 = 1, on the other hand if A2 7^ 1, twinned martensite laminates appear. However, we did not find significant differences 
between the transformation hysteresis of Ti 50 Ni27Pd23 and Ti 50 Ni39Pdii alloys. Indeed, as it is explained above, the behavior of 
hysteresis in thermoelastic martensites depends, in principle, on frictional mechanisms. It is experimentally observed that when 
A2 = 1, the austenite/martensite interfaces are essentially coherent and display very few dislocations ifTTl . This is probably due 
to the fact that, when A2 = 1, the compatibility between austenite and martensite may be fulfilled with a simple habit plane, i.e., 
without any transition layers associated to unrelaxed elastic energy. Generally, when A2 7^ 1, incompatibility enforces complex 
habit planes between polytwinned martensitic domains and the austenite, and it may be that the associated stored elastic energy 
is relaxed by plastic accommodation along the habit planes. Recently, this behavior was experimentally observed in reference 
[6]. As a result, it may be expected that internal friction of the moving interfaces is significantly higher when \ 2 ^ 1, resulting 
in a hysteresis larger than in the case A2 = 1. 

In conclusion, the reason why we could not observe the differences in behavior of hysteresis of different alloys is understood 
by the fact that our phase-field model does not incorporate the dynamical dislocation creation during the transformation. The 
relation of the width of the transformation hysteresis to the value of A2 requires a modeling that couple the phase transformation 
with plastic activity (58l[59l- Finally, our results suggest that the behavior of hysteresis width cannot directly be understood by 
assuming that the condition A2 = 1 decreases only the critical nucleus size (see [ 38 ]). It is probably a combined effect together 
with the decrease of friction at the interfaces when \ 2 = 1. This conclusion has also been supported by recent numerical studies 
(39j [40j, where the generation and accumulation of persistent lattice defects during forward-and-reverse phase transformations 
are observed. 
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Appendix A 

The elastic energy of a coherent inclusion of a single variant of martensite (variant p) embedded into the austenite, can be 



calculated using Eq. 1 1 



lasti' 



Ijf B w (n)Mq)| 2 ^L (Al) 



The matrices B pp are calculated by inserting the deformation strain tensor (p) of variant p into Eq. 12 Since the term B pp (n) 
is positive, when the order parameter r\ p is different than 0, the minimum of the strain energy is obtained for a platelet having a 
volume V aligned in the dir ectio n of n that minimizes the function B pp (n). We can numerically calculate the direction no that 



minimizes the strain energy ( Al ). Calculations show that it exists a habit plane normal n for which the matrices B pp (n) vanish 
for the alloys Ti 5 oNi3gPdn and TiNiCu satisfying the condition A2 = 1. On the other hand, the matrices B pp (n) ^ for any 
value of n, when X 2 ^ l.lt means that austenite/monovariant habit planes do not cost any strain energy only when \ 2 = 1. 

In order to calculate the elastic energy of a coherent inclusion of a poly-twinned domain, formed by two variants p and q, 
embedded into the austenite, we calculate the average deformation of the poly-twinned martensite domain 

€ij = K€ij(p) + (1 - K)€ij(q), (A2) 



where k is the twin ratio. We can substitute the average deformation into Eq. [AT] and minimize it with respect to k, to find 
the twin ratio and habit plane direction no of the twinned-martensite domain with the austenite. We found that, when \ 2 = 1 
condition is fulfilled, the elastic energy approaches zero when twin ratio k goes to zero. 
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